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Abstract. In a Randall-Sundrum type II inflationary scenario we compute 
perturbation amplitudes and spectral indices up to next-to-lowest order in the slow- 
roll parameters, starting from the well-known lowest-order result for a de Sitter brane. 
Using two different prescriptions for the tensor amplitude, we show that the braneworld 
consistency equations are not degenerate with respect to the standard relations and we 
explore their observational consequences. It is then shown that, while the degeneracy 
between high- and low-energy regimes can come from suitable values of the cosmological 
observables, exact functional matching between consistency expressions is plausibly 
discarded. This result is then extended to the Gauss-Bonnet case. 
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1. Introduction 

The advent of new Kaluza-Klein models whose features can be tested by modern 
astrophysical, cosmological and ground-based collider observations has triggered interest 
in those scenarios in which our visible universe is embedded in a higher-dimensional 
spacetime, with extra dimensions that can be 'mildly' compact or even non-compact 
[I] See, e.g., [HE] for a review. In particular, the Randall-Sundrum type II model 
[6j has received much attention since its birth [H] PB] because of its rich theoretical 
structure in a relatively simple conceptual framework. One of the many aspects of 
the model is the possibility for the cosmological evolution to change substantially 
from the standard four-dimensional case, since the extra non-compact dimension can 
communicate via gravitational interaction with the matter confined in the brane 
universe. Gravity is free to propagate in the bulk, which is assumed to be empty in 
the simplest scenarios. 

According to the present data, the large-scale structure of the universe, as well as 
the anisotropies of the cosmic microwave background, can be explained by an early stage 
of accelerated expansion driven by a slow-rolling scalar field, which boosted primordial 
quantum fluctuations up to cosmological magnitude. The consequent perturbation 
spectra have been widely studied in the standard 4D scenario thanks to a number 
of techniques, including the slow-roll (SR) formalism we will adopt here (for a general 
review and notation, see [13); i n particular, it is possible to combine the available 
observables in relations that are characteristic of the inflationary paradigm and verify 
them through CMB and sky-survey data. 

Braneworld calculations for the perturbation spectra are much more involved 
because of the complicated geometrical tissue and only quite general formalisms or 
approximated approaches have been explored so far. In this paper we continue the 
discussion of the inflationary consistency equations started in j!5j . in which next-to- 
lowest SR order expressions for perturbation amplitudes and indices are constructed 
from the lowest-order, de Sitter results known in the literature. Here we will give 
another prescription for the tensor amplitude motivated by recent 5D results [16J, 
showing however that there are no new observational features with respect to the original 
ansatz. The reliability of the consistency equations for testing high-energy effects is then 
addressed; it is shown that experimental degeneracy can come from a low tensor-to-scalar 
amplitude ratio as well as a nearly scale-invariant spectrum. However, simple slow-roll 
considerations involving the SR structure of the cosmological large-scale observables 
show that exact functional degeneracy of at least one relevant consistency equation 
is expected to be discarded, both in RS and Gauss-Bonnet gravity, whatever the full 
five-dimensional amplitudes may be. 

The plan of the paper is the following: section El sets up the background model; 
sections |3] and |U are devoted to the SR treatment of the cosmological perturbations; in 
section El and in the appendix we study the consistency equations and their structure, 
while the final discussion is in section |U1 
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2. Basic setup 

The basic model consists in a Randall- Sundrum type II scenario with a brane- confined 
inflaton field and no energy exchange between the brane and the bulk. The continuity 
equation is the standard 4D equation coming from the local conservation of the stress- 
energy tensor, V y T^ u = 0: 

p + 3#(p + p) = 0, (1) 

and the modified Friedmann equations in this conformally flat background are |9~] [TB"] 

H 2 = ^P(2X + p) , (2) 

- = -^[A(p + 3p) + p(2p + 3p)], (3) 
a da 

where H is the Hubble parameter restricted to the 3-brane, k\ = 8ir/m1 includes 
the effective Planck mass WI4 ~ 10 19 GeV, A is the brane tension (constrained to 
be A 1 / 4 > 100 GeV by gravity experiments) and p is the matter density, which is 
p = cf> 2 /2 + V = p + 2V in the case of a homogeneous scalar field with potential 
V{4>). The case of brane-bulk energy flow is far richer and will not be treated here (e.g., 
[T7]-j2n]). In the following, superscripts (I) and (h) will denote the low- (p A) and 
high- (p ^> A) energy regimes, respectively. 
The slow-roll parameters are 

(4) 
(5) 

^ V 2 ■ (6) 




H 2 ct) 

Since a/a — H 2 (l — e), the condition for inflation is precisely e < 1 and inflation ends 
when e = 1. The time derivatives of the SR parameters are 

e=2He(^e- v y (7) 

77 = H (£77 -£ 2 ), (8) 
where = 1 and q^ = 2. For future reference we note that, in the high-energy limit, 

<P 2 = ^He= P -e, (9) 
where p = a/ Ak|/6 . 

3. Scalar perturbations 



The 5D Einstein equations for a brane with an isotropic fluid embedded in an anti-de 
Sitter bulk are very complicated due to both the great number of degrees of freedom for 
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the cosmological perturbations and the non-local physics coming from the possibility 
for KK gravitational modes to propagate and interact throughout the whole spacetime. 
For this reason, while a general setup has been established [211 |301 ? a ^ uu 5D spectrum 
amplitude has not been calculated yet for either scalar or tensor modes, except in the 
case of some particular scenarios, for example those in which the brane is de Sitter 
[121 EI] [US] or in the large-scale limit [2B 123 I2H EH]- In this context the scalar 
perturbation is found to be the same as the one in the standard four-dimensional 
background. A way to see this in the slow-roll framework is to note that in the large- 
wavelength limit, k <C aH, the four- dimensional perturbed metric induced on the brane 
can be expressed by means of the two longitudinal-gauge-invariant scalar potentials $4 
and \l/4, 



ds 



brane 



1 + 2$ 4 ) dr - a{ty{l - 2^ 4 )% dx'dx 3 . (10) 



In the case of brane-bulk interaction through a non-diagonal stress-energy bulk tensor, 
the complete continuity equation is p + 3H(p + p) + 2tb = 0, where tb is the 05- 
component of the bulk tensor (see [23] )• Then the 5D equations of motion for both the 
homogeneous background scalar field <po{t) and its fluctuation 50(t,x) are 

0o + 3#0 O + V'ifa) + 2^ = , (11) 

00 

5<f> - a- 2 V 2 5(f) + 3H5(f) + V"5<f> - ($5 + 3* 5 )0o + 2$ 5 V + 4$ 5 ^ = , (12) 

00 

where $5 and ^5 are the gauge-invariant scalar potentials which appear in the full 5D 
metric and a prime denotes a 0-derivative. By setting an empty bulk (r# = 0) one 
recovers the ordinary 4D equations and in the large-scale limit $5 — ^5 ~ = $4 — ^4. 
At this point one can conjecture that the perturbed effective gravity and scalar 
actions giving the equations of motion are the usual 4D ones, and performs pure four- 
dimensional calculations, for example by the methods of [HZ] . Since the explicit form of 
H is not involved, as one can easily check, one obtains the standard Mukhanov equation 



z" 



<+ (^ 2 --J Mk = 0, (13) 
where primes denote second derivates with respect to conformal time 

r = /7 = -^(M' (14) 

and «k are the coefficients of the plane wave expansion of the canonical variable 
u = —zlZ. Here, z = a<p/H and 1Z is the curvature perturbation on comoving 
hypersurfaces. Generally, the spectral amplitudes are defined as 

257H z z 

in particular, the scalar amplitude reads [3D| 

2 v-3/2 T(u) / l H 



A* 



5tt T(3/2) 



k=aH 
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where v = a/1/4 + t 2 z" / z ~ 3/2 + 0(e) is a combination of SR parameters and the 
expression is evaluated at the horizon crossing of the perturbation with wavenumber k. 
This amplitude comes from an exact solution of the Mukhanov equation with constant 
SR parameters. A group of braneworld solutions found in jH] will suit well since they 
approach the power-law inflation [32] at late times ^H] (see also |43j). Perturbing around 
one of these solutions and assuming that e, r\ and £ are small enough (see equations (J7|) 
and (jHJ)), it follows that 



S 1 ' 2 = A S = [l-(2tf + l)e + CW 

071 



1 H 



(16) 

k=aH 

where C ~ —0.73 is a numerical constant; by fixing the square bracket equal to 1 one 
gets the lowest-order expression. As a final remark we note that this is not a true 5D 
derivation, nor has the full five- dimensional backreaction of the metric been taken into 
account. So equation (fT7)|) must be regarded as an effective result rather than an exact 
one. Moreover, for simplicity we have neglected any contribution coming from both the 
total anisotropic stress and the Weyl tensor (and their fluctuations); however, we would 
expect this extra-physics not to dramatically modify the above inflationary perturbation 
spectrum, since it would only affect the small-scale/late-time structures (HHl E31 EEl HEj ■ 
The spectral index is defined as 

din S 1 S , irA 

^- 1 "dinT = F(i^)5' (17) 

where we have used the exact relation dln/c/d0 = H(l — e)/0. To first SR order, the 
scalar index is uniquely defined in both regimes, 

n 5 - 1 ~2r?-4e, (18) 

while the energy dependence is absorbed in the second-order part via equation [THj . 
Defining the running as = dn^/dln/c, it is easy to see that 

af =2(-4e 2 + 5er / -a, (19) 
af = 2(-2e 2 + 5e V - £ 2 ) . (20) 

In order to write these expressions in terms of observables, we will assume that \rj\ 
or, alternatively, that £ does not change in the energy scale, ~ 



4. Tensor perturbations 



In the case of a de Sitter brane (p +p — 0), it is possible to decouple the zero-graviton 
mode from the massive tower of Kaluza-Klein gravitons, thus obtaining a first-order SR 
expression for the tensor amplitude J3*T] : 

2 



T 1 / 2 = An 



^f(* 



(21) 



k=aH 



where 



Fix) 



Vl + x 2 - x 2 \n ll/x + a/1 + 1/x' 



-1/2 
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In the limit p <C A, F(x) ~ 1 and one recovers the standard spectrum, while in the limit 
p 3> A, F(x) ~ a/3s/2. One way to state this result is that, to lowest-order SR, the 
braneworld tensor amplitude is given by the 4D expression under the mapping 

ho: H i-> HF(H/fi) . (22) 

The amplitude (|21|) and its scalar counterpart give rise to a consistency equation which 
is the same as the standard 4D relation; in order to check whether the degeneracy is an 
accidental feature of this particular approach, higher-order amplitudes and consistency 
equations are needed. Since an enhanced calculation in a non-decoupled regime is a very 
difficult task, truly 5D-computed expressions have been performed only in a variety of 
limits and approximations, none of which to second order SR. 

In [16J a perturbed de Sitter brane whose Hubble constant is experiencing a 
discontinuous variation SH = Hi—H 2 is studied. A 5D spectral amplitude T 5D , obtained 
as an expansion in 5H/H <C 1, is provided at mildly large scales, k <C H 2 /5H, and is 
compared with the 4D braneworld amplitude T e ff derived from a generalization of the 
map ()22j) acting on the standard amplitude T^o oc H 2 , namely, 

hi: HF(H/p) . (23) 

Here H = H[1+0(5H/H)] includes the correction to the lowest-order result coming from 
the 4D zero-mode-zero-mode Bogoliubov coefficients. It turns out that the lowest-order 
0(1) and next-to-lowest order 0(SH/H) effective amplitudes match the five-dimensional 
result, 

rp(0) — rp(0) _ rp(0) rp(l) _ rp(l) (r,A\ 

1 — 1 eff — 1 5D J 1 eff ~ 1 5D 5 

and \ (T 5 D — T" e //) /T^°'\ ~ 0(5H J ' H) 2 outside the horizon, while inside the horizon there 
can be significant discrepancies. 

The case of a continuous smooth variation of the Hubble parameter, which is typical 
during the inflationary regime long before the reheating, has to be treated separately 
because of the time dependence of H in the equations of motion of the KK modes. 
However, one expects the slow-roll generalization of the previous result to display 
similar features ^H], since e = 0{5H / H). In fact, in order the variation SH not to 
be damped by the accelerated expansion, it must occur in a time interval 8t ~ H^ 1 , 
that is, H ~ —5H/5t ~ —5HH, giving the heuristic correspondence e ~ 5H/H. Thus 
we can conjecture an analogous mapping (f2~3J) in the SR expansion such that holds, 
with H = H[l + 0(6,77)]. We schematically represent it in figure d 

Here fo maps any next-to-lowest-order function to its lowest-order SR formt. We 
remark once again that, while f Q is trivial and go = 1 was first demonstrated in [31], it 
would be a non-trivial goal to show that g\ = 1 for a non-de Sitter inflationary brane 
or, more generally, to find the region in the space of the effective amplitudes where /o 
is one-to-one. 

f This is equivalent to setting e = r\ = only in the amplitude T, because, for instance, fo(S) oc tT x . 
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T, 
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Figure 1. The /i-mappings lll'l'l) and <\'2:>li between 4D and effective 5D tensor 
amplitudes. 



There is however another possibility. According to equation (fTH|) . the mapping (|2^|) 
can be regarded as acting on z rather than on the Hubble parameter: 

h: z i — ► z/F(H/n) , (25) 

where z = z^d — \/2a/ k± for the tensor perturbation. This more elegant formulation has 
two advantages: first, it encodes all the needed SR information in one single relation, 
so h\ and h,2 collapse to each other; secondly, it is reasonably valid in more general 
braneworld scenarios [TFj . 

In , the following ansatz for the tensor amplitude to next-to-lowest order in SR 
parameters was proposed: 



A T =-^[l-(C+l)e]— f(- 



(26) 



k=aH 

It is straightforward to see that both the low-energy regime [33] and the extreme SR 
regime [31] are regained. From the above discussion it emerges that another prescription 
is viable for the (squared) Randall-Sundrum tensor amplitude: 

(27) 



T = T 4D f*°Y 



k=aH 



\ Z J 

In the high-energy limit, 

T W= [l-(3C + 2)e](3 2 H 3 \ k=aH , (28) 



where 



02 _ 1 



25tt 2 4/i ' 

and we have used the approximated function F(x) ~ \j3x/2 together with equations 
(J7J), (fTlj) and (j2Zj). The tensor index reads 

n ^ = H(i- 7) W) = ~ e[3 + {3C + 5)e ~ 2{3C + 2)r]] ' (29) 

to be compared with the low-energy standard result 

n§? = -2e[l + (2C + 3)e-2(C + l)77]. (30) 



Degeneracy of consistency equations in braneworld inflation 



8 



5. Consistency equations and braneworld degeneracy 



5.1. Experimental degeneracy 

In the lowest-order approximation, the only available consistency equation is degenerate, 



that is, n 



(h) 



n 



(I) 



-2T/S 113113. Using the relations 



[l-2C(e-77)]| 



2T 



for p <C A 
for p ^> A 



e ~ [1 - C(e 

we find from equations JTBJ), (fTHjl - ^ . (|27|l . (J2T3J) and (|5Dj| . 



n 



CO 



2 T - 

S 

T 



' S 



l-| + (l-n 5 ) 



2T 2, 

1-33 + 3(1-".) 



a 



(0 



a 



T 

2T 

3S 



L2- + 5(n 5 -l) 



32 T 



(31) 
(32) 

(33) 
(34) 
(35) 
(36) 



Equation (J5l"j) confirms the idea that the new prescription is more correct than that of 
[T5] . since it bears no trace of the computational method; conversely, in equation (54) 
of [12] the constant C, typical of the solution of the Mukhanov equation, does appear 
explicitly. 

We do not write the consistency equation for the running of the tensor index since 
it is degenerate and therefore does not give any additional information. By combining 
the last four equations, one gets 



Aut 
Aa s 



Tl/rp 71/rp 



(h) 



a 



— a 



(h) 



2T 
3S 
IT 

3S 



I + (n S ~ 1) 
'44 T 



(37) 
(38) 



Note that, provided T / S + (ns — 1) >0 and as < 0, the high-energy spectrum has 
a redder tensor tilt and a steeper scalar running. It was pointed out in j^lj that the 
failing of the standard lowest-order consistency equation is due to the introduction of the 
formal limit A — >• 00 , which decouples any bulk physics from pure 3-brane quantities but 
does not necessarily lead to a "smooth" connection between 5D gravitational quantities 
and their 4D equivalent. Indeed the tensor amplitude is a typical example of this 
mechanism, its new features being encoded both in the explicit if-dependence and in 
the SR structure. 

In order to test the observational signatures that the braneworld scenario may 
produce, we will use the estimates of the cosmological observables given by the first- 
year analysis of the WMAP mission [SHj^ISS] (see also the BOOMERanG experiment 
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|q"4] and the Ly-a forest observations [55]). Substituting the upper bound for the tensor- 
to-scalar amplitude ratio T/S = 0.06 jJU] and the best fit for the scalar index of [52J, 
ns = 0.93, we find from the previous equations 



which is close to the error in the estimate of jHSl, &s — —0.04 ± 0.03, as well as to 
the accuracy estimate of the future Planck mission [SB]. The effect coming from the 
consistency relation for the tensor index is too weak to be detected by the missions of 
this generation: putting ny — —0.1, equation (pTTjl gives —Aut/ut = O(10~ 4 ), which 
does not improve the result coming from equation (|26|) [To] . 

Of course the most practical way to obtain experimental degeneracy is to consider 
functionally different consistency equations, whatever they are, in the two regimes 
and small tensor-to-scalar ratios as well as a nearly invariant spectrum (ns — 1), 
an eventuality which is quite possible in the range of the currently available data. 
For instance, taking the ratio T/S = (T / S) max /2 = 0.03, which is within the 2a 
likelihood bound of |SZj, we have Acts — 0.004 and —Aut/ut = 0.006, while for 
T/S = [T/S) max /3 = 0.02 we get Aa s ~ 0.002 and -An T /n T = 0.003, an order of 
magnitude smaller than the most optimistic high-ratio case. Figure El shows a three- 
dimensional plot of equation (|3*H|) as a function of the amplitude ratio and the scalar 
index, while in figure 01 some specific values of ns are considered. Note that the break of 
the degeneracy increases for higher values of the spectral index; moreover, the relative 
sign of the effect may change when ns < 1. 



Aa s ~ 0.01, 



(39) 




i 



0.08 0.9 



Figure 2. Aas as a function of the ratio T/S and the scalar index ns- 




Figure 3. Aas for particular values of ns- From bottom to top (increasing thickness), 
n s = 0.9, 0.95, 1, 1.05 and 1.1. 



5.2. Theoretical degeneracy 

At this point one may wonder if there can exist next-to-leading-SR-order perturbation 
amplitudes that give the same consistency equations for the standard four- dimensional 
cosmology, thus ruling out the possibility to discriminate between the two scenarios. 
However, to construct mathematical expressions without any physical content would 
be of little use since they would not provide a physical explanation why consistency 
equations should be still degenerate; so the consistency equation approach would still 
be worth investigation. After this preamble, it should be noted that it is not possible to 
construct simple perturbation amplitudes that give degenerate consistency equations 
and reproduce the extreme SR limit (j21j) . The reason is the following. We will 
concentrate on the consistency equation involving the scalar running as since it is the 
most relevant relation from an observational point of view. It is an equation coming from 
the lowest- SR-order part of the scalar index, equation (fT%|) . which is generated by the 
lowest-order part of the scalar amplitude, that is the functional part of (fTBj) . We recall 
that each time derivative raises the SR order by one at any step: the time derivative of 
the functional part of the amplitude gives the linear part of the spectral index, while 
another derivation gives the running, through equations (J7J) and (jHJ)- Therefore one 
should change the functional part of the amplitude S in order to impose = ctg\ 
regardless of the explicit SR structure. Since equation (fTBj) has not been computed 
in a truly 5D framework and the bulk source term of (jllj) has not been taken into 
consideration, this could be feasible a priori, even if equation (frfljl is well supported in 
many respects^. An example of a generalized scalar amplitude is (the superscript (h) 

| We stress once again that the lowest-order scalar amplitude can be calculated, in the long wavelength 
part of the spectrum, by a number of different background- independent techniques. 
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S=[l + f s {e, V )]— sl -, (40) 



and evaluation at horizon crossing are understood) 

1 25?r 2 t 

where b and c are constants, s is a normalization pre-factor and fs is a linear function 
of the SR parameters and of any dimensionless combination of cosmological quantities 
such as m 4 , A, H and its time derivatives. 

As has been said in the previous section, the mapping h is well motivated only in 
the quasi-de Sitter case, so in principle a more general high-energy structure for the 
tensor amplitude than that of equation (}2Tj) is possible, keeping the functional part 
constrained by the zero-mode 5D calculation: 

T=[l + f T (e }V )](3 2 H\ (41) 

where fx is the tensor counterpart of fs- For reasons of simplicity, we have dropped 
the H dependence in the SR functions since it gives rise to a polynomial structure for 
S which does not change the main argument; see the appendix for details. Also, for the 
above considerations we can ignore these SR functions. 
Equation (j4*Uj) generates the spectral index 

n s — 1 ~ br] — ce (42) 

and the running 



as ^ e[— ce + (b + 2c)rj\ = e 



2c\ „ , 2c 



2 



(43) 



The case b = ^ c (c = ^ b) is discarded because it is not possible to absorb the rj (e) 
dependence of the running into a cosmological observable. The standard case c = 2b = 4 
is trivial. Now, assuming 
T 

- ~ re , (44) 

r being a constant, equation (|33|) holds in the high-energy limit if b = 24r 2 /(5r — l) 2 
and c = 12r 2 /(5r - 1) . Using equations 0, ©, (HH and flUJ) we have, to lowest SR 
order, 

T ~ e !#4-^&-2 e b/2 p b/2-c + 3 (45) 

S 2s 



this implies that a high-energy relation such that equation (|44|) is possible only in the 
standard case. In general, a traditional consistency relation will not be obtained in a 
braneworld context and the brane tension will appear explicitly, as well as the inflaton 
potential. Suitable powers of these quantities will be set by imposing that the ratio (pToj) 
is dimensionless: if d s is the dimension of s, then it must be d s + c — 2b = 0. 
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5.3. Generalization to patch cosmology and the Gauss-Bonnet case 

We can generalize the previous arguments to patch cosmology and, in particular, to the 
Gauss-Bonnet (GB) scenario. In this case the tensor index consistency relation is no 
longer degenerate (see and jSHj for details and references). The modified Friedmann 
equation of the patch approach is H 2 oc p q , q being the same constant as in equation 
(J7J), and <f) 2 oc pe, as before. For the GB high-energy regime, q = 2/3. Equation 
becomes 



1 + (n s -l)+(e-l + ^)ce 



(46) 



where 6 = 2(1- q' 1 ) = -1, 1, for q = 2/3 (GB), 2 (RS) and 1 (4D), respectively; the 
tensor amplitude is T oc H 2+e . One may also try to compensate by hand for the e factor 
of (|45|) with a slightly more general scalar amplitude. Looking at the ratio (|45|). it is 
clear that there could be a scalar amplitude able to absorb the e and p factors, namely, 

S = — -e-^ , (47) 
25n 2 6 b 



where 7 is a constant. The spectral index and its running are given by equations (|42|) 
and (J4T)|) when substituting b and c with 

b -> b + 27 , 

c — > c + 7 . 

In particular, equation (|43|) now reads 

l OC e 7+V2 p [2(2 9 -l) +6 - gC ]/2 . (4g) 

in order to satisfy equations (|33j) and (jSj), we must solve the system 

b i 
7 + ^ = 1, 

2(2q- 1) + b-qc = 0, 

c + 7 
6 + 2 7 

(#-2 + 5r) (c + 7) = 12r 2 . 

For general q, there do exist non-trivial solutions. In the Gauss-Bonnet case, however, 
r = 1 in equation f)44j) . as in four dimensions [3*H] . and again the system admits no 
solution; the same considerations apply for the Randall-Sundrum model, r = 3/2, 
regardless of the introduction of the new degree of freedom 7. 



6. Discussion 



A new era of high-precision cosmological observations (WMAP, 2dF [59J, SDSS 
|6T1] and others) can now constrain the high-energy physics of braneworld models 
predicting significant deviations from the standard four-dimensional big bang scenario 
|3H1 EB1 E3 EI]^|E11- In this perspective, it is useful to investigate what are the 
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main features that can discriminate these models from the 4D case; in particular, an 
early inflationary period experienced by a Randall-Sundrum brane can display different 
spectra of cosmological perturbations. 

In this paper we have investigated the relations linking the cosmological observables 
generated by this accelerated expansion. By extending the "/i-mapping" of the de 
Sitter case to a braneworld scenario with a continuously and smoothly varying Hubble 
parameter, it is possible to construct a reasonable next-to-lowest-order expression for the 
tensor amplitude starting from the 4D amplitude of the same order; this procedure leads 
to a set of consistency equations which carry informations about the extra dimension 
and that can be probed by the cosmological experiments of this generation. It turns out 
that, even if the SR structure of the observables is rather rigid, the consistency equations 
of the standard four- dimensional inflation no longer hold in the brane cosmology and 
the lowest-order SR degeneracy is softly broken by 5D gravitational effects. 

Therefore, while observational degeneracy of the consistency relations is achievable 
within the range of cosmological parameters determined by recent experiments, the 
theoretical structure is unstable against even long-wavelength 5D contributions (via the 
effective Friedmann equations (j2J) and (JHJ)). With explicit examples, we have shown that 
it would be easy to spoil the simple closed structure of these equations while driving 
off the standard scalar and tensor high-energy amplitudes, both for Randall-Sundrum 
and Gauss-Bonnet scenarios. More accurate SR calculations could show consistency 
relations with a different structure, perhaps unveiling terms that cannot be directly 
related to the available observables [31] • Whatever the final answer turns out to be, 
there is hope of tracking down braneworld signatures through large-scale inflationary 
physics. 
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Appendix: slow-roll structure of cosmological observables 

In this section we provide a general argument to show that the lowest-order linear 
structure of the spectral indices and the quadratic one for their running are a 
consequence of the slow-roll approach, thus justifying the use of a scalar amplitude 
S oc H c /(f) b as a meaningful probe in analysis. This, however, does not say anything 
new about the consistency relations because their closed structure depends on the form 
of the tensor-to-scalar ratio, say, equation flHj) . First of all, from equations (jlj) and (JZJ) 
we have 




(A.l) 



Now define a tower of slow- roll (hereafter, 'nSR') parameters as 

d n+i H 

e " = ~Hh>H ' 



(A.2) 
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where d n H indicates the nth time derivative of H. In general, this definition would not 
be convenient for two reasons. The first is that, when expanding around a flat potential, 
the convergence of the parameters to zero would be made unclear by the derivative in the 
denominator. The second reason is that it is not possible to truncate any expansion to a 
finite slow-roll order because all the nSR, parameters are involved with the same power, 
unlike what happens with the standard definitions. However, since we are interested 
in the behaviour of leading-order expressions in the sense of equations (0}-©, for our 
purpose (|A.2j) will be sufficient. See [HH] for a systematic study of the SR tower in the 
four- dimensional case. Equations (@J, © and (|A.1|) tell us that 

e =e , (A.3) 
2 V = -0eo + ei, (A.4) 

where 9 = 2(1 — q^ 1 ), while, noting that 

e n = He n (e + e n - e n+l ) , (A. 5) 

H = ffVi , (A.6) 

from (|SJ) we have 

2£ 2 = eie 2 -e? + #eo(eo-ei). (A.7) 

In this nSR formalism, the explicit energy dependence is removed from the complete 
slow-roll tower. 

In general, a cosmological amplitude A can be written as a sum of dimensionless 
combinations of H and its derivatives: 

A = J2 Ak > ( A - 8 ) 

k 

A k = a k \{ (d»»H)- l »» , (A.9) 

n k £U k 

where a k is a normalization pre-factor, U k C N and the exponents l nk £ 14 C Q are 
determined by the physics of the inflationary period. The spectral index is defined as 
ha = (ha.A)'/H, that is, 

n A = E * Ak f , (A.10) 

l^k A k 

(fc) _ 1 A k 

n\ ' = . A. 11 

A HA k 1 ; 

Using the distributive property for generic I?jS, (n^i) = Y^i &i Ilj^ > after some 
steps one gets 

n k n k 

According to this formula, by the very definition (jA.lOJ) the spectral index is linear in 
the nSR parameters e; when l Uk ^ for some n k , since 0(ua) = O(n^) . By equations 



(k) 
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i"3|) and (jA.4j) . this means that the lowest-order part of the spectral index is at least 
linear in e and 77 provided that there are non-zero coefficients§. 
The running of the spectral index is, using ()A.10|) and (jA.llj) . 

aA= H-^AT + ^AT- n ^ (A ' 13) 



where 



(fc) _ n A 



l ^n k {e + e„ fc - e„, fe+ i) . (A.14) 

Again, we find a rigid slow-roll structure that constraints the running to be 0(n\) or, 
in the standard SR parameters, to be at least 0(n A ) to leading order. 
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